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ABSTRACT 



An infinite integral over four spherical Bessel functions of the form 

r"^ j\^{kir) j\^{k2r) jx^{kir) j\^{k2r) dr is analytically evaluated. The result 
involves a finite sum over an associated Legendre function of integer degree and half 
integer order with a real argument greater than 1. Evaluation of such functions is 
discussed. 
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1. Introduction 



Infinite integrals over spherical Bessel functions have always been of interest due to 
their occurence in nuclear physics [1-10], particle physics [11] and astrophysics [12-13], 
to mention a few. In this paper, an infinite integral over four spherical Bessel functions 
of the form 



is analytically evaluated for the special case when = ki and k^ = ^2- The 
parameters Ai, A2, A3 and A4 are assumed positive integers. Recently, a generalised form 
of the integral has been attempted numerically [14], and analytically [15]. However, 
the analytical evaluation involves a complicated hypergeometric function. In our 
evaluation, we start by assuming the momenta ki, /c2, ^3 and A;4 to be positive and 
form the sides of a quadrilateral, which is a consequence of the conservation of linear 
momentum in any physical application. The general result of evaluating eq. (1.1) is 
in the form of a finite sum over 3j and 6j symbols with a finite integral over Legendre 
functions with a rational function. This integral is evaluated for the special case when 
ks=ki and k4=k2 in appendix A. The final result for this special case involves a finite 
sum over an associated Legendre function of integer degree, half-integer order and a 
real argument greater than 1. Appendix B discusses the evaluation of such functions. 
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2. Evaluating The Infinite Integral Over Four Spherical Bessel Functions 
The integral 

oo 

(2.1) 



where ki, k2, ks and k4 are positive real numbers, can be written as 
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X yj r''^ jxsihr') jx^ikAv') jLiKr') drj , 
using the spherical Bessel functions Closure Relation 

oo 

J K^jL{Kr)jL{Kr')dK = ^5{r' r). (2.3) 



The value of L is chosen such that it is the smallest value that satisfies 

|Ai - A2I < L < Ai + A2, (2.4) 

and 

IA3 - A4I < L < A3 + A4. (2.5) 
Previously [7-10] it was shown that 
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where A = (kj + k^ - K'^)/2kik2 and ^(A) = ^(1 - A)^(l + A) with 9 the Heaviside 



function in half- maximum convention. Pi{x) is a Legendre polynomial of degree I, 

Ai A2 L\ fAi A2 l1 

is a 3i symbol and < > is a 6i symbol which can be 

oy \c L-c I ] 

found in any standard angular momentum text [17-18]. Similarly 
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where A' = {kj + kj - K^)/2ksk4. The result for eq. (2.1) is then 
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, , > J{ki,k2,k3,k4;l,l ,L). 



00 

J{kuk2MM;l,l',L) = j ^i^}^^Pi(A)Pi,{A')dK. (2.9) 



This integral can be analytically evaluated for the special case when ks = ki and 
ki = k2, i.e. A' = A as follows: 
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Transforming the variable of integration in (2.9) from to A the integral becomes 
.ft, fe, ... kr, I, I', L) = / .A, (2.10) 



where y = {kf + k2)/2kik2. 

Using appendix A, (2.10) reduces to 



J{ki,k2,ki,k2;l,l',L) ^" '^'^^^ 



r(L + l/2) \kl - kl\L 



(2.11) 



The associated Legendre function, P^^ ( |fc^-fc^| )' integer degree and a 

half-integer order with a real argument that is greater than 1. Appendix B discusses 
the associated Legendre function for such cases and shows closed form expressions for 
special cases of L. Furthermore, if A2 = Ai = li, A4 = A3 = I2, the ideal choice for L 
is L = 0. Eq. (2.1) then becomes 

Jr''jUkir)ji,{k2r)jUkir)jUk2r)dr = ^Y.(^^^ (2.12) 
where A;< (/c>) is the smaller (larger) of ki and /c2- 
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3. Conclusions 



The integral over four spherical Bessel functions 

oo 

(3.1) 



was analytically evaluated for positive integer Ai, A2, A3 and A4. The resulting 
expression involved a finite sum over an associated Legendre function with integer 
degree and half-integer order. An interesting compact result was found upon setting 
A2 = Ai and A4 = A3. 
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APPENDIX A: Evaluation of the Integral ^^^d^dA 
Starting with the integral from [19], eq. 7.228, page 830 

L— 1/2 

where (y) is a Legendre function of the second kind of degree /j, and order L — 1/2. 

Also, using [19], eq. 8.739, page 1023 gives 

-i^(L-i/2) 0^-1/2/ X _ /^ r(L + ^ + 1/2) ^^_i/2 y 
Hence, the itegral in {A.l) becomes 

1 

f Pf^jx) _ /g- r(L + // + l/2) 2 ^^/2 p-M-1/2. y _ ^ 3. 

Multiplying eq. (A. 3) by (2// + 1) P^(A) and summing over /j, results in 



(^.4) 



using 

^(2// + l)P^(A)P^(x) = 2S{A- x). (A5) 

Hence 
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-p^=) y (2L' + 1) ^ ^^Pi,(A), 



{A.6) 



using 



pm(A^pM(.^_(_-,^m+M I {f^ + m)] {,y + M)] 
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{A.7) 

So, multiplying eq. {A.6) by Pi'{A) and integrating over A from —1 to 1 results in 



j fa_A).-.l/2 ''^ = '^ r(L+l/2) EP/'+i) 
-1 



(A.8) 
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APPENDIX B: Associated Legendre Fumctions, Pp{x) for a; > 1 

Previously we have shown that associated Legendre functions, P^'^"(a;), for integer I 
and real m can be written as 

for real —l<x<l, and the factorial is replaced by the gamma function to allow for 
non-integer m. The special cases for particular I values were 

= ^ TV (15x='-15x2m-9a: + 6xm2+4m-m')(ii^)'»/^ (B.5) 

J- ( o 7n\ ~\~ 1) 1 3y 

Pn^) = — N (105 x"^ - 105 x^m - 90 x"^ + 45 xV + 55 xm 

r( 4 — m + 1) ^ ^ 

- lOxm^ + 9 - 10 + m^) (1±:^) 

1 — X 

Now, by comparing equations 8.702, page 1014 and 8.704, page 1015 of [19], one finds 
that for real x > 1 

with the special cases 
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^."'w= r(|i -Li + (I^)"^^- f^-') 
= r(i2-H+i) (3x^-3x™-l+,.')(i±|^/^ (B.io) 

^3"(a;) = T^THt; — ^ {lbx^-lbx^m-9x + 6xm^+4m-m^){^^)"'/'^, (S.ll) 

r(|3 — m\ + 1) X — 1 

P^{x) = — — — — (105a;^- 105 x^m- 90x^ + 45 a; V + 55 



lOxm^ + 9 - 10 + m^) (:i^)W2^ 

a; — 1 
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